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Abstract: 

Investigations of neutron(n)-capture element nucleosynthesis and chemical evolution have 
largely been based on stellar spectroscopy. However, the recent detection of these elements 
in several planetary nebulae (PNe) indicates that nebular spectroscopy is a promising new 
tool for such studies. In PNe, n-capture element abundance determinations reveal details of 
.s-process nucleosynthesis and convective mixing in evolved low-mass stars, as well as the 
chemical evolution of elements that cannot be detected in stellar spectra. Only one or two 
ions of a given trans-iron element can typically be detected in individual nebulae. Elemental 
abundance determinations thus require corrections for the abundances of unobserved ions. 
Such corrections rely on the availability of atomic data for processes that control the 
ionization equilibrium of nebulae (e.g., photoionization cross sections and rate coefficients 
for various recombination processes). Until recently, these data were unknown for virtually 
all )i-capture element ions. For the first five ions of Se, Kr, and Xe — the three most widely 
detected rc-capture elements in PNe — we are calculating photoionization cross sections 
and radiative and dielectronic recombination rate coefficients using the multi-configuration 
Breit-Pauli atomic structure code AUTOSTRUCTURE. Charge transfer rate coefficients are 
being determined with a multichannel Landau-Zener code. To calibrate these calculations, 
we have measured absolute photoionization cross sections of Se and Xe ions at the Advanced 
Light Source synchrotron radiation facility. These atomic data can be incorporated into 
photoionization codes, which we will use to derive ionization corrections (hence abundances) 
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for Se, Kr, and Xe in ionized nebulae. Using Monte Carlo simulations, we will investigate 
the effects of atomic data uncertainties on the derived abundances, illuminating the systems 
and atomic processes that require further analysis. These results are critical for honing 
nebular spectroscopy into a more effective tool for investigating the production and chemical 
evolution of trans-iron elements in the Universe. 

PACS Nos.: 32.80.Fb, 33.60.+q, 34.80.Lx, 34.70.+e, 95.30.Dr, 95.30.Ky, 97.10.Cv, 97.10.Tk, 
98.38.Bn, 98.38.Ly, 98.38.Hv 

1. Introduction 

1.1. Astrophysical Motivation 

The formation of neutron(n)-capture elements (atomic number Z > 30) is not as well understood 
as that of lighter species, and observational abundance determinations are critical for testing and im- 
proving theories for the nucleosynthesis of these elements. Trans-iron elements are predominantly 
formed via rapid or slow n-capture nucleosynthesis (the "r-process" and ".s-process," respectively), in 
which the distinction is based on the relative timescales of n captures and (5 decays. Each mechanism is 
responsible for the production of roughly half of the n-capture element isotopes in the Universe [1,2]. 

Trans-iron elements are sensitive probes of the nucleosynthetic histories of astrophysical objects 
and the chemical evolution of the Universe, in spite of their low abundances [2]. Much of our knowl- 
edge of the synthesis of these elements is based on the interpretation of abundances derived from stellar 
spectra [1-3]. However, stellar spectroscopy reveals only part of the story of how n-capture elements 
are produced, due to the limited number of elements accessible in stellar spectra, and the difficulty of 
observing stellar photospheres during certain stages of evolution. In recent years, nebular spectroscopy 
has been shown to be a potentially effective tool for investigations of n-capture element nucleosynthe- 
sis, based on the detection of these species in a large number of planetary nebulae (PNe 4 ) [4-7]. 

The progenitor stars of PNe (1-8 M Q ) may experience .s-process nucleosynthesis during the late 
asymptotic giant branch (AGB) phase of evolution, which precedes the formative stages of PNe. In this 
process, called the "main" .s-process to differentiate it from the "weak" .s-process that operates in mas- 
sive stars [8, 9], free neutrons are produced primarily by the reaction 13 C(a, n) le O in the "intershell" 
region between the H- and He-burning shells. Fe-peak nuclei experience a series of n captures and 
/^-decays to transform into heavier elements. In more massive AGB stars (> 4-5 M Q ), the intershell 
region may reach temperatures allowing the 22 Ne(a, n) 25 Mg reaction to efficiently produce neutrons. 
The intershell material, enriched in .s-process nuclei and carbon, can be convectively dredged up to the 
stellar atmosphere and ejected into the interstellar medium by stellar winds and PN ejection [1, 10]. 

Nebular spectroscopy provides unique information regarding the nucleosynthesis of n-capture ele- 
ments, opening a window to elements that are not detectable in cool giant stars. For example, nebular 
spectroscopy provides access to Ge, Se, Br, Kr, Rb, Te, I, Xe, Ba, and Pb, each of which has been iden- 
tified in the spectrum of at least one PN [4—7, 11-15]. With the exceptions of Rb, Ba, and Pb, these 
elements are not detectable in AGB stars, and in fact had not been detected in any of their sites of origin 
until their identification in PN spectra. Therefore little is known about the origins and chemical evo- 
lution of these elements beyond theoretical predictions that have not been tested against observational 
measurements. 5 Moreover, nebular spectroscopy enables n-capture nucleosynthesis to be investigated 
in classes of stars and stages of evolution during which the stellar photosphere is not easily studied. 
Intermediate-mass (^4-8 M Q ) AGB and post- AGB stars experience high mass-loss rates that enshroud 
their atmospheres in optically thick, dusty circumstellar envelopes [19, 20]. As a consequence, little 



4 Following astrophysical nomenclature, we abbreviate the singular form planetary nebula as "PN" and the plural 
form as "PNe." 

5 However, recent models of AGB nucleosynthesis have been applied to Se and Kr abundances in PNe [16-18]. 
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is known about their contribution to the Galactic inventory of n-capture elements. However, these 
objects are observable as Type I PNe (characterized by large N and He enrichments as expected for 
intermediate-mass stars [21, 22]), in which n-capture elements have been detected [6, 17]. Similarly, 
once an AGB star becomes carbon-rich (C/O > 1) the large increase in opacity leads to heavy mass 
loss [23, 24], obscuring the photosphere with circumstellar material. Convective dredge-up during the 
final thermal pulses present computational difficulties for AGB star models, but can have an important 
effect on n-capture element yields [17]. The ionized portion of PNe are composed of material enriched 
by the final thermal pulses during the AGB, and hence provide an excellent venue for constraining the 
dredge-up efficiency at low envelope mass. Finally, n-capture elements are potentially excellent tracers 
of C enrichments in PNe, as both are formed in the intershell region of AGB stars and are convectively 
dredged up simultaneously. Given the difficulty in accurately determining C abundances in ionized 
nebulae [25, 26], the s-process enrichments of trans-iron elements in PNe provide sensitive probes for 
enrichments of this astrophysically and biologically important element. 

1.2. Atomic Data Needs 

Before the full potential of nebular spectroscopy of trans-iron elements can be realized, several 
important atomic data problems must be addressed so that n-capture element abundances can be accu- 
rately determined in ionized nebulae. One must first determine the abundances of the observed ions, 
and then correct for those of unobserved ions (wherein the greatest uncertainties lie). Much of the 
atomic data needed for these determinations are unknown, due to the recent nature of this field (n- 
capture elements were not identified in any ionized nebula until 1994 [4], and only in the last few years 
have they been detected in a large number of nebulae - totaling about 100 objects to date [5, 6]). 

For ionic abundances, the atomic data needed include wavelengths, transition probabilities, and 
effective collision strengths. For most of the n-capture element ions detected in nebulae, energy levels 
are satisfactorily known, though the NIST compilations [27] are less complete than for lighter elements. 
One notable exception is Se 2+ , which can have significant fractional abundances particularly in low- 
ionization nebulae. The fine-structure energies of the Se 2+ 3 P ground term are uncertain by 6.3 cm -1 
[28], which lead to a wavelength uncertainty of -10 A for the transition [Se III] 1 D 2 - 3 Pi 8854.2. 
Notably, this line was reported to be detected in the PN NGC 7027 [4], and Sharpee et al. [5] noted 
that the nearby He I 8854. 14 line showed marked excess relative to other He I lines in the same series 
in two PNe. More accurate energies for the Se 2+ 3 P fine-structure levels are clearly needed. 

Forbidden line transition probabilities have been determined for the ground configurations of most 
detected n-capture element ions. These were computed by Biemont & Hansen for ions with 4p n 
valence shells [29-32] and by Biemont et al. for 5p n valence-shell ions [33]. In contrast, effective 
collision strengths for electron-impact excitation have been calculated only for Kr 2+ -Kr 4+ [34] ; Xe 2+ , 
Xe 3+ , Xe 5+ , Ba+, andBa 3+ [35]; and Se 3+ andZr 3 " 1 ". 6 To our knowledge, effective collision strengths 
have not been computed for other n-capture element ions, and are particularly needed for Br and Rb 
ions. 

In a typical ionized nebula, only one or two ions of a given n-capture element can be detected 
(due in part to the low cosmic abundances of these species, but also to the electronic structure of each 
ion and the accessibility of strong emission lines). The conversion of ionic abundances to elemental 
abundances therefore requires ionization correction factors (ICFs) for the abundances of unobserved 
ions. The most robust method of determining the relevant ICFs is through numerical simulations of a 
nebula's thermal and ionization structure. However, such models rely on the availability of accurate 
data for atomic processes affecting the ionization equilibrium of these elements. For photoionized 
nebulae such as PNe, these data include photoionization (PI) cross sections and rate coefficients for 



6 K. Butler, unpublished 
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radiative recombination (RR), (low-temperature) dielectronic recombination (DR), and charge transfer 
(CT) reactions. 7 Until recent work by our group, these data were not known for the vast majority 
of trans-iron element ions, thereby rendering n-capture element abundance determinations in ionized 
nebulae uncertain by an estimated factor of two or more [14]. 

In the remainder of this paper, we concentrate on recent work that will provide the atomic data 
needed to determine ICFs and hence elemental abundances of Se, Kr, and Xe - the three most widely 
detected «-capture elements in PNe and other ionized nebulae. We have utilized theoretical methods to 
calculate PI cross sections and rate coefficients for RR, DR, and CT rate coefficients, and benchmarked 
our results with experimental absolute PI cross section measurements. We discuss our theoretical work 
in §2 and the experimental measurements in §3, before providing concluding remarks in §4. 

2. Theoretical Atomic Data Results 

We have computed PI cross sections and RR, DR, and CT rate coefficients for ions of Se and Kr, 
and similar calculations for Xe are in progress. While these data are in principle needed for each of the 
> 30 ions of each element, in practice we are able to simplify the problem to the neutrals and first five 
ions, since PN central stars are not hot enough (generally < 2 x 10 5 K [36, 37]) to significantly ionize 
species with ionization potentials greater than 100 eV. 

For PI, RR, and DR, we have used the atomic structure code AUTOSTRUCTURE [38], which can 
efficiently compute the electronic structure, radiative and autoionization rates, and multi-configuration 
distorted-wave PI cross sections for each ion. Relativistic effects are accounted for via Breit-Pauli 
formalism and semi-relativistic radial wavefunctions. RR rate coefficients were determined from the 
direct PI cross sections using detailed balance, and DR rate coefficients from the resonant portion of 
the PI cross sections. The independent processes approximation [39] was used to treat RR and DR 
separately. The results presented in this section are preliminary, and will be fully detailed elsewhere. 

The electronic structure was computed with Thomas-Fermi-Dirac-Amaldi model potentials, and 
the average of the LS term energies were optimized by varying the radial scaling parameters. The con- 
figuration interaction (CI) expansions and term energy optimizations were adjusted to best reproduce 
the experimental energy levels and ionization potentials cited in NIST [27]. 

We computed PI cross sections out of the valence shell for all levels in the ground configuration 
of each ion, using the CI expansions from the structure calculations. Ground state PI cross sections 
for Se and Kr ions near their ionization thresholds are presented in Figure 1. RR rate coefficients 
(total and partial final-state resolved) were then computed over the temperature range (10 1 -10 7 )2: 2 K, 
where z is the charge. For RR, the CI expansion of the N+l -electron ion was expanded to include 
all one-electron additions to the target configurations. Total and partial final-state resolved DR rate 
coefficients were computed over the same temperature range. For this process, we focused on An = 
core excitations, since high-temperature DR is expected to be negligible at photoionized temperatures. 
For DR, the target CI expansions were augmented in order to account for all relevant core excitations, 
and 1 -electron additions to these configurations were included in the N+ 1 -electron ion CI expansions. 
Following Badnell [40], we used experimentally-determined target energies whenever possible for our 
DR calculations, using the theoretical level splittings in cases where the NIST energies are incomplete. 

We estimated the uncertainties of our results, as these can have significant effects on nebular abun- 
dance determinations. We used three CI expansions of different sizes for each ion, and also tested the 
sensitivity of our results to other internal parameters of the code (e.g., scaling parameter values for 
continuum orbitals, the use of non-orthogonal orbitals, the internal interpolation scheme, etc.). We 
also compared our results to experimental PI cross sections (see §3). 



7 At photoionized nebula temperatures, processes such as collisional ionization and high-temperature DR negligi- 
bly affect the ionization equilibrium. 
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Fig. 1. Calculated PI cross sections for the ground state of each Se and Kr ion considered, in the energy region 
near the ground state ionization threshold. 



For PI, we estimate typical uncertainties of 30-50% in the direct cross sections near the ground 
state ionization threshold. The largest uncertainties come from using different CI expansions, and are 
largest for low-charge species. Similar uncertainties are found for our RR rate coefficients. 

Uncertainties in low-T DR rate coefficients primarily arise from the poorly known energies of 
low-lying autoionizing levels, as only these states are available as DR channels at the temperatures 
of photoionized nebulae. The near-threshold resonance positions are notoriously difficult to calculate 
precisely, leading to large uncertainties in the low-T DR rate coefficients. Unfortunately, the energies 
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of these near-threshold states have not been experimentally determined for n-capture element ions (or, 
indeed, for most row 3 and heavier elements [41, 42]), and hence it is unknown whether they lie 
just above the threshold (and are available as DR channels) or just below, in which case they do not 
contribute to DR (though our use of observed target energies alleviates this effect to some extent). 8 To 
test the sensitivity of our rate coefficients to this effect, we shifted the continuum relative to the near- 
threshold resonances by an amount equal to the largest energy difference in theoretical vs. experimental 
energies. These tests revealed uncertainties in our rate coefficients near 10 4 K (the typical electron 
temperature of PNe) ranging from a factor of 2 to an order of magnitude for most Se and Kr ions. 

At 10 4 K the DR rate coefficient is greater than that of RR by 1-2 orders of magnitude for low- 
charge Se and Kr ions, and hence low-T DR is the dominant recombination mechanism. In photoion- 
ization codes used to model ionized nebulae, the large uncertainties in low-T DR rate coefficients 
necessitates that this process be treated approximately. As an example (by no means an egregious 
one!), Cloudy [44] assumes that the low-T DR rate coefficients for the first four ions of Row 3 and 
heavier elements are the same as the average rate coefficients for C, N, and O ions of the same charge. 
However, given the complexity of the target states for near-neutral n-capture element ions, there is a 
much denser set of near-threshold resonances that can contribute to low-T DR than for L-shell sys- 
tems. Thus, as atomic number increases the approximation in Cloudy becomes increasingly poor, and 
for Se and Kr the computed low-T DR rate coefficients exceed the averaged CNO rate coefficients by 
1-2 orders of magnitude. The accuracy of low-T DR rate coefficients is a challenge common to all 
photoionization codes, and is important to our understanding of the ionization balance of astrophysical 
nebulae. 

CT requires a different treatment, as it is a quasi-molecular problem. We consider collisions with 
H atoms of the form X q+ + H — > x( 9_1 )+ + H+ + SE, where q = 2-4, since this process is not 
expected to play a major role in the ionization equilibrium of more highly-charged ions due to the 
minimal mixing of neutral H with those species in PNe. We are using a multichannel Landau-Zener 
code based on the formalism of Butler & Dalgarno [45] and Janev et al. [46] to determine total and 
partial final-state resolved CT rate coefficients for Se, Kr, and Xe atoms. This method has been shown 
to be accurate to within a factor of ^3 for lighter elements [45, 47]. Singly-charged ions, however, 
cannot be treated in this approximation as there are no avoided crossings in the adiabatic potentials 
(CT is a Demkov process for those species) [48]. For these systems, we will use a modification of the 
approximate procedure of Pequignot & Aldrovandi [49] to estimate the CT rate coefficient. 

3. Experimental Results 

As is true of any atomic data calculation, it is important to benchmark the results to experimental 
measurements. To this end, we have measured absolute PI cross sections of Se and Xe ions at the 
Advanced Light Source (ALS) synchrotron radiation facility at the Lawrence Berkeley National Lab- 
oratory in California. PI cross sections for Kr [50-52] and Xe 3+ -Xe 6+ [53, 54] had previously been 
measured, and hence we focused our efforts on the first five Se ions and the first two of Xe. 

These measurements were performed using the merged-beams technique [55] at the Ion Photon 
Beamline (IPB) apparatus [56, 57] located on beamline 10.0.1 of the ALS. In these measurements, 
a plasma of atomic ions is created within an electron-cyclotron-resonance source and is accelerated 
with a potential (typically 6 kV). The ion of interest is selected with an analyzing magnet to form the 
ion beam, which is collimated with two sets of vertical and horizontal slits and focused with three 
electrostatic einzel lenses. A set of spherical electrostatic bending plates merges the ion beam with 
the photon beam, and along the merged beam path photoionization occurs. Photoions are redirected 



Recent work [43] suggests that the situation may not be as dire as described here. In finite density plasmas, the 
continuum threshold is broadened due to thermal considerations, and resonances just below the canonical ioniza- 
tion threshold can contribute to DR. Thus, the associated uncertainties in the rate coefficient may be decreased. 
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to a detector by a demerging magnet and counted. The photons are produced by electrons in the 
1.9 GeV storage ring of the ALS that are accelerated by a 10-cm period undulator. Photon energies 
are selected with a spherical monochromator with three adjustable gratings to cover energies between 
18 and 300 eV. The energy resolution is controlled by adjusting the entrance and exit slits of the 
monochromator. 

By measuring the photoion yield as a function of photon energy, the photoionization spectrum is 
determined in arbitrary units. To place this on an absolute scale, the photoions must be produced in a 
well-defined volume. This is achieved by placing a potential on the "interaction region," a stainless- 
steel mesh cylinder with a precisely determined length. The potential energy-tags photoions produced 
in that region, and the demerging magnet setting is altered to direct only energy-tagged photoions to 
the detector. The volume in which the photoions are produced is determined by measuring the overlap 
of the photon and ion beams with three translating-slit scanners in the interaction region. Absolute 
measurements are performed at discrete photon energies to normalize the photoionization spectrum. 

We have measured absolute PI cross sections near the ground state ionization thresholds of Se + 
(Figure 2, left panel), 9 Se 2+ , Se 3+ (Figure 2, right panel) 10 , Se 5+ , Xe + , and Xe 2+ . The Se 4+ cross 
section could not be measured due to a large background also seen for the isoelectronic Kr 6+ (caused 
by the presence of ions in autoionizing metastable states in the primary ion beam). Measurement 
and identification of the resonance structure near the ground state threshold is valuable, since these 
resonances correspond to low-T DR resonances and can be used to constrain DR rate coefficient cal- 
culations. 

Most atomic ions in photoionized nebulae are in the ground state, and thus PI out of the ground state 
is of most interest for astrophysical applications. However, the absolute PI measurements at the ALS 
are complicated by the fact that the primary ion beam is not composed solely of ground state ions, but 
contains significant fractions of metastable states with lifetimes exceeding the time of flight through the 
IPB apparatus. The metastable fractions can be constrained by measuring the PI cross section below 
the ground state threshold down to the threshold of the highest energy metastable state (for example, 
the ground state ionization threshold of Se + is 21.2 eV; in Figure 2 the cross section and resonances 
below that energy are all from metastable states), and with the aid of theoretical calculations. 

4. Concluding Remarks 

We have presented preliminary results of an investigation to determine atomic data needed for 
accurate abundance determinations of the three most widely-detected n-capture elements in ionized 
astrophysical nebulae - Se, Kr, and Xe. We are utilizing the atomic structure code AUTOSTRUCTURE 
[38] to calculate PI cross sections and RR and DR rate coefficients for the first five ions of these 
elements, and are using a MCLZ code to determine CT rate coefficients. To calibrate our theoretical 
results, we have measured the absolute PI cross sections of Se and Xe ions at the ALS synchrotron 
radiation facility. We also make use of existing Kr and Xe absolute PI measurements [50-54]. 

The end goal of determining these new atomic data is to utilize them to derive analytical ionization 
correction factors (ICFs) for the abundances of unobserved Se, Kr, and Xe ions in ionized nebulae. 
Following the methodology of [14], we plan to incorporate the new atomic data into the photoioniza- 
tion code Cloudy [44], which we will use to run a grid of models over a range of stellar and nebular 
parameters (e.g., stellar effective temperature and luminosity, and nebular density) typically encoun- 
tered in PNe. ICFs can be derived by finding correlations between the fractional abundances of the 
observed Se, Kr, and Xe ions and those of routinely detected species such as O or Ar ions, and fitting 
the correlations with analytical functions. The new ICFs will be applied to existing spectra of n-capture 



'Results have been submitted to J. Phys. B: At. Mol. Opt. Phys. for publication, 
"festeves et al., in preparation 



©2010 NRC Canada 



8 



Can. J. Phys. Vol. XX, 2010 




£ 5 



o 

O 2 
1 





18 19 20 21 22 23 24 25 26 27 
Photon energy (eV) 



42 44 46 48 50 52 54 56 
Photon energy (eV) 



Fig. 2. Left panel: The Se + absolute photoionization cross section measured at the ALS. The solid line is the 
photoionization spectrum normalized to absolute cross-section measurements at discrete photon energies (circles 
with error bars corresponding to the 90% confidence level). Resonances and Rydberg series are indicated. Right 
panel: The Se 3+ absolute photoionization cross section, with absolute measurements indicated by circles with 
error bars indicating the 90% confidence level. This figure is modified from [7]. 



elements in PNe [5-7, 14] to reassess the abundances of these three elements. This will permit robust 
abundance determinations of n-capture elements in PNe (and other types of photoionized nebulae), 
enabling investigations of s-process nucleosynthesis in their progenitor stars at an unprecedented ac- 
curacy. 

We will also test the sensitivity of abundance determinations to our estimated atomic data uncer- 
tainties using Monte Carlo simulations. This will help to assess which systems and atomic processes 
require further analysis. Additionally, this will highlight the general dependence of abundance deter- 
minations on the underlying atomic data that forms the foundation of astrophysical spectroscopy. 
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